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Turbulence modeling 

By T. -H. SHIH 

1. Motivation and Objective 
Motivation 

In recent years codes that use the Navier-Stokes equations to compute aerody- 
namic flows have evolved from computing two-dimensional flows around simple 
airfoils to computing flows around full scale aircraft configurations. Most flows 
of engineering interest are turbulent and turbulence models are needed for their 
prediction. Yet, we know that present turbulence models are adequate only 
for simple flows and do poorly in complicated flows such as three-dimensional 
separation, large-scale unsteadiness, etc. The same progress that allowed the 
development of these aerodynamic codes, namely the introduction of supercom- 
puters, has allowed us to compute directly turbulent flows, albeit only for simple 
flows at moderate Reynolds numbers. These direct turbulence simulations pro- 
vide us with detailed data that experimentalists have not been able to measure. 
This work is motivated by the fact that data exists for developing better turbu- 
lence models and by the need for better models to compute flows of engineering 
interest. 

Objective 

The objective of this work is then to develop turbulence models for engineering 
applications. The model categories that show promise for immediate use are on 
the two-equation level and the Reynolds-stress level. We will make use of existing 
methodologies to develop models. The models will be tested using data from 
direct simulations, experiments and analysis. Specifically, our objectives are as 
follows: 

1. Examine the Reynolds stress budgets using direct simulation flow fields 
(Mansour et al 1988, Moin et al 1989). 

2. Use Rapid Distortion Theory to analytically study the effects of mean 
deformation on turbulence. In particular, examine the development of the rapid 
pressure-strain under rapid distortions. 

3. Compare existing models with data and theory. Develop models where 
needed using appropriate expansions and constraints. Test these new models 
using results from direct simulation, experiment and theory. 

4. Use the method of moment generating function to extend second order 
closures to higher order closures. We know that there exist a close connection 
between high-order moments and coherent structures. The moment generating 
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function approach should allow us to include the effect of coherent structures on 
the models. 

5. Extend turbulence modeling to compressible flows. 


2. Work Accomplished 

1. Numerical simulation of a three-dimensional boundary layer — AIAA pa- 
per no. 89-0373 (work supported in part by AFOSR) 

The objective of this simulation is to study the mechanics of three-dimensional 
boundary layers and develop improved models for their predictions. Three- 
dimensional effects were achieved by direct simulation of a fully developed tur- 
bulent channel flow subjected to transverse pressure gradient. To obtain a good 
statistical sample during the transient period, 14 computer runs were ensem- 
ble averaged. Each run started from a different realization of the channel flow 
(far apart in time). The simulation shows that, in agreement with experimental 
observations, the Reynolds stresses are reduced and that near the wall a lag 
develops between the stress and the strain rate. The reduction in the stress is 
due largely to a drop in the production rate and an increase in the dissipation 
rate. In the coming year, we will study the performance of existing second order 
closures in predicting this drop and introduce improvements where needed. 

2. k-e modeling 

k-e turbulence models are often used in computing engineering flows with 
moderate success. In general, the flow field is qualitatively well predicted, but 
quantitative agreement often falls short. In the case of Reynolds stress modeling 
the question of a length scale for an eddy viscosity does not arise, and the energy 
dissipation rate, e, is one of the unknowns of the problem. In this work an e 
equation has been developed following the methodology of Lumley (1978). We 
used it in connection with both two-equation modeling and Reynolds stress 
modeling. 

A two-equation model for use in low-Reynolds number flows has been tested 
against the channel data of Kim, Moin and Moser (1987) . The modeled equations 
are given by, 
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o Jfc = 1 
— 1.3 
Ci= 1.48 
C 2 = 1.8 
e = € — (}> 
tj* = kjkj/ (2fc) 
i*r = Cpfpk 2 /*. 

C„= 0.09 

/„ = 1 - exp(-0.0115— ) 



We find that the model adequately predicts the mean velocity profile (see figure 
l) and the turbulent kinetic energy (see figure 2). 

3. Second order modeling of near-wall low-Revnolds number turbulence 

A set of second order closure models for low-Reynolds number turbulence 
has been developed for the simulation of wall bounded flows without using wall 
functions. The wall effect is built in the pressure-strain correlation term of the 
Reynolds stress equation and in the modeled terms of the dissipation rate equa- 
tion. We find that realizability is particularly important for modeling the near 
wall turbulence. The proposed models are particularly suitable for surfaces of 
arbitrary topology since they do not use the wall distance as a parameter. The 
models are tested by computing the fully developed channel flow. The full set of 
equations are used to compute the mean velocity, all the Reynolds stresses and 
the dissipation rate of the turbulent kinetic energy. We find reasonable agree- 
ment between the prediction of the mean profile and the data (see Figures 3) 
which is an indication that the shear stress is well predicted. However, the nor- 
mal stresses are not as well predicted. In particular the streamwise component 
is underpredicted (see figure 4) while the transverse component is overpredicted. 
The cause of this shortcoming is still being investigated. 

4. Second order modeling of a passive scalar in turbulent shear flows — AIAA 
paper no. 89-0607 

A model equation for the scalar dissipation rate was proposed using the ansatz 
that the ratio of mechanical time scale to scalar time scale has an equilibrium 
value. In addition a model for the pressure related terms in the scalar flux equa- 
tion was constructed based on consideration of realizability. The models were 
tested by comparison with experimental data for heated plane and axisymmetric 
jets. 

5. Rapid Distortion Theory and 2-D 2-C turbulence modeling 

Rapid distortion theory was used to analyze the development of the Reynolds 
stress in a 2-D 2-C homogeneous turbulence under mean irrotational strain and 
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mean shear. Here, 2-C refers to two-component turbulence. We found that the 
development of the Reynolds stress under 2-D 2C conditions is very different 
from its development under 3-D 3-C conditions. The findings can be summarized 
as follows: 

1. The mean shear or strain have no effects on the isotropy of the Reynolds 
stress if the initial field is isotropic. 

2. Initially anisotropic fields become isotropic under the influence of mean 
shear or strain. 

These findings are opposite to the finding in 3D 3C turbulence where mean 
shear or strain will drive isotropic turbulence away from isotropy. In 2D 2C 
turbulence the pressure strain drives the Reynolds stresses back to isotropy. 

A general 2D 2C model for the pressure-strain term, Tii, was constructed 
which contains only one undetermined coefficient C, 

2^ = Si, / "H ^CbjkbpjSpic + 2fl jkbik + 2n ik b k ] 

where, 

Si,- = \ (U,.i + Vi.,), n,,- = i {V,J - Vi.,) 

, u7u7 1 _ , 1 

bij -~n- 3 %. k ~ 2 U ’ U ’ 

Rapid distortion theory will be used to determine this coefficient. We postulate 
that the more general 3D 3C model should reduce to this model in the limit of 
2D 2C turbulence. 

3. Future Plans 

1. Examine the performance of existing second order closures in predicting 
three-dimensional turbulent boundary layers. Introduce improvements where 
needed. 

2. Use Rapid Distortion Theory to study the effects of mean deformation on 
turbulence. Extend the analysis developed for mean strain and shear to study 
the effect of rapid rotation on the turbulence. 

3. Use the method of moment generating function to extend second order 
closures to higher order closures. 

4. Use non-weighted ensemble averaging method (Shih et al, 1987) to extend 
the second order models to compressible flows. This averaging technique (as op- 
posed to Favre averaging) will allow us to apply all the incompressible modeling 
methodologies to modeling compressible flows. 
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Figure 1. Mean velocity profile in a channel as predicted by a two-equation 
model. 













